The SED in the hot continuum of the 
symbiotic binary ARPavonis 



I. Tests with current models 



cn 
o 
o 

(N 
X) 

CN ! Abstract 



> 



A. Skopal^ 



Astronomical Institute, Slovak Academy of Sciences, 059 60 Tatranskd Lomnica, 

Slovakia 



We present the spectral energy distribution (SED) in the continuum of the echpsing 



fSJ \ symbiotic binary ARPav between 0.12 and 3.4 /um. This revealed a high luminosity 

■ of the hot object in the binary, ~ 2 200 ((i/4.9kpc)^ Lq. We introduce a method 
. of disentangling the total continuum spectrum into its individual components of 

^) \ radiation for current models of symbiotic binaries. Applying a standard ionization 

■ model we show that the configuration of ARPav differs significantly from that 
typical for symbiotic binaries during their quiescent phases. The best fit of the 



observed SED is provided by radiation of a simple blackbody accretion disk with 



Lad ~ 1 700 ((i/4.9kpc)^L0, which is embedded in an extended hot corona with 
2 ■ Te = 40 000 ±5 000 K and Ln ~ 500 (ci/4.9kpc)2L0. This basic configuration of 

0^ \ the hot object explains also the observed wavelength-dependent depth and width 

of the eclipse profile. The standard thin disk model requires a high accretion rate 
-^acc ^ 2 X 10~^ Mq yr~^ onto the central star with a radius -Race ^ 2 Rq to balance 
^ , the observed luminosity. Even irrespectively to the disk model, the accretion process 

^ ; limits i?acc > 0.1 i?0 for Mace > 1-7 X W-^ MQyr-^ and Mace = 0.75 - 1.0 Mq, 

which precludes a white dwarf to be the accreting star. Application of models with 
the disk and the boundary layer shows that the far-UV spectrum is not consistent 
with a large amount of a hot radiation from the boundary layer. However, the 
presence of such a boundary layer in the system is indirectly indicated through the 
strong nebular emission. This solution suggests that the hotter inner parts of the 
disk including the boundary layer are occulted by the disk material in the direction 
to the observer. 

Key words: Eclipsing binaries, Symbiotic stars. Accretion and accretion disks 
PACS: 97.80.Hn, 97.80.Gm, 97.10.Gz 



E-mail: skopal@ta3.sk 
Preprint submitted to New Astronomy 2 February 2008 



1 Introduction 



AR Pavonis was discovered W iMavalj (Il937h as an eclipsing P-Cygni type star 
with a period of 605 days. AndrewsT l 1974 ) found that the eclipsed object is 
highly v ariable in both brightn ess and size. On the basis of ultraviolet spec- 
trosco py Hutchings et al.l (|l983[ ) confirmed this view. Bruch. Niehues and JonesI 
f)l994 ) observed a periodic (~ Porb) wave-like variation, which occurred in the 
visual light curve between eclipses from 1985.7. As the most satisfying expla- 
nation they suggested a brightness modulation due to variation in the mass 
transfer rate from the red giant to t he hot component. Analysing historical 
1889-1998 light curve, ISkopal et al.l (|2000ar i identified the real change of the 
orbital period at a rate P = —3.5 ± 0.8 x 10~^ between 1896 April and 1985 
September. They also noted that during quiescence, the hght cu rve has char 



acteristic features similar to those observed in CVs. Recently, ISkopal et al. 
(j2001a|) refined the mean orbital period to 604.45±0.02 days. 



As to the ARPav configuration, Thackerav and Hutching^ ( 1974) proposed 
a binary model, in which the red giant fills its Roche lobe and loses mass 
to an evolved, ~30 000K hot, compact object surrou n ded b y a moderately 
dense gas cloud or a thick ring. iKenvon and WebbinkI ()l984 ) found that the 
observed UV continuum colours for ARPav are well comparable to those de- 
rived fror n a model o f an ac cretion disk around a main sequence accretor. 
Recently, ISchild et al.l (|2001l l. having new radial velocities of the red giant, 
determined its mass function to 0.055 and the projected rotation velocity 
of ll±2kms~^. Assuming co-rotation, they derived the giant radius of 130 7?© 
and by comparing its adequate stellar parameters with RGB/AGB evolution- 
ary tracs in the HR diagram, they suggested the giant mass as 2Mq, and, 
consequently from the mas s function, the hot component mass as 0.75 Mq. 
Based on these parameters, ISchild et al.l (|200lf ) suggested that the hot com- 
ponent in ARPav is most probably a white dwarf and the cool component is 
deeply inside of its Roche lobe (however, as we will show later in this paper, 
the giant in th i s syst em does most probably fill its Roche lobe). Most recently, 
Quiroga et al1(|2002l ). having radial velocities of both components, determined 
their masses as Mg = 2.5 Mq, Mh = IMq. They obtained the radial velocity 
curve of the hot star under assumption that the broad emission wings of the 
Ha, H/3 and Hen 4686 profiles originate from the inner accretion disk or an 
envelope around the central star. In addition, the authors suggested that the 
central absorption in hydrogen line profiles is formed in the neutral part of 
the giant's wind, which is concentrated in the orbital plane. 



In this paper we analyse the SED in the continuum between 0.12 and 3.4 yum 
with the aim to determine a basic configuration of the hot object in ARPav. 
In Sect. 3 we introduce individual components of radiation often used in inter- 
pretation of the spectra of symbiotic stars. In Sect. 4 we compare the observed 
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SED and those given by current models of symbiotic stars, including a stan- 
dard ionization model and models with an accretion disk. Finally, in Sect. 5 
we discuss their physical plausibility and suggest basic structure of the hot 
object. 



2 Observations 



To reconstruct the SED of ARPav in the range of 0.12 to 3.4 /im we used 
the broad-band optical and infrared UBVRI and JHKL photometry and the 
lUE (International Ultraviolet Explorer) low-resolution spectroscopy. 



The optical/IR magnitudes were summarized bv ISkopal et al.l fcOOOal) and 
some new observations were recently published by Sko pal et al. (|2000b[ ). To 
de termine t h e star 's bri ghtness in eclipses we used photometric measurements 
of Andrewsl ()l974l ) and lMenzies et al.1 (|l982h . Observations in the i?c and /c 
bands of the Cousins system were transformed into the Johnson system accord- 
ma; to iBesselll (|l983h Stellar magnitudes were convert ed to fluxes according 
to the calibration of iLena. Lebrun and Mignar dl (|l999l). 



The ultraviolet fluxes were taken from Fig. 5 of ISchild et al. ( 2001 ). We used 
average values from those measured out of eclipses (i.e. between orbital phases 
0.91 > ip > 1.1; the rms within 10-15%) and values from the mid eclipse. To 
get bette r coverage of the ultrav iolet continuum we added a few values from 
Fig. 2 of iHutchings et all (Il983h and measurements presented in Table 4 of 
Kenvon and WebbinkI (|l98J). 



Finally, we dereddened observations with -E-r-v = -26 (|Kenvon and Webbink 
1984 ) using the extinction curve of iMathisI ( 1987 ) in the UV/optical region. 



while in the IR wavelengths we used reddening ratios of Savage and Mathia 
f)l979>) . Table 1 summarizes the observations used in our paper. 



3 Spectral energy distribution 



Figure 1 shows the 0.12 - 3.4 /im SED in the continuum of ARPav. This 
suggests a high luminosity of the hot object, similar to that of its cool com- 
panion. The out-of-eclipse continuum is characterized by: (i) A round far-UV 
shape with decreasing fluxes to shorter wavelengths, and (ii) a fiat profile in 
the near- UV/optical region. Both the steep far-UV (A < 1500 A) continuum 
and the prominent Balmer jump in emission - characteristic features of the 
continuum of symbiotic stars during quiescence - are not present. We will fit 
the out-of-eclipse continuum by that including an accretion disk, but also test 
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Fig. 1. Dereddened fluxes of the AR Pav continuum given by the low-resolution lUE 
spectra, optical UBVRI and infrared JHKL photometry (Table 1). Measurements 
made during eclipses are marked by (*), out of eclipses by (•), a nd open circles 



(0) c orrespond to fluxes of the giant derived from its spectral type (|Menzies et al 

Observations during eclipses can be approximated by nebular continuum 
of Te = 40 000 K (N) and fluxes of the giant by a Planck function of the colour 
temperature Tg = 2 990K (G). The dashed thick line represents their sum. 

whether a standard ionization model of symbiotic binaries can explain these 
characteristics. 



During eclipses, the ultraviolet continuum weakened by a factor of ~3-5, it 
is flatter than that observed out of eclipses and displays a discontinuity at 
U/B. This suggests the nebular-type spectrum of a very hot gas. We will fit 
the uneclipsed continuum by the nebular radiation of a high electron temper- 
ature. In the following sections we describe simple models of the considered 
components of radiation, which we use to fit the observed SED. 



3.1 Nebular continuum 



In the conditions of ionized gaseous nebulae in symbiotic stars (a dense medium 
with Ue ~ 10^ — 10^^ cm~^), the main contributor to the near-UV/optical con- 
tinuum is the free-bound (f-b) emission from hydrogen. Less important are 
the f-b emission of He ii and He i. Due to the dominance of the H i f-b emis- 
sion and its basically simila r profile with that of He i and He ii (see Fig. 1 of 
Brown and Mathewsl . 1970l ). and also for the sake of simplicity, we calculate 



here only Hi f-b and f-f nebular contributions. The nebular flux is proportional 
to / Uj^Uq dV, where n+ and Ue is the concentration of protons and electrons, 
respectively. Accordingly, the observed flux, F^^^ (ergcm~^ s^^ A^^), of the 
nebular continuum at the wavelength A, can be written as 



A-Kd^Ff"' = ex J n+ne dV, 
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Table 1 

Continuum fluxes between 0.12 and 3.4 ^m. 



A[A] 


-^obs 




-^obs 




Ref. 




in eclipses 


out of eclipses 




1270 


0.15 


1.36 


0.40 


3.75 


3 


1300 






0.48 


4.27 


4 


1350 


0.23 


1.85 






3 


1380 


0.25 


1.87 


0.64 


4.67 


5 


1455 


0.26 


1.83 






3 


1520 


0.25 


1.70 


0.73 


4.96 


5 


1600 


0.27 


1.76 






3 


1700 






0.91 


5.88 


4 


1808 


0.28 


1.80 


1.07 


6.87 


5 


1875 


0.26 


1.67 






3 


2200 


0.20 


1.94 


0.50 


4.88 


3,4 


2400 


0.23 


1.30 






3 


2600 


0.29 


1.40 


0.80 


3.89 


3,4 


2700 


0.25 


1.13 


0.93 


4.22 


5 


3000 


0.34 


1.28 


1.59 


6.00 


5 


3100 


0.40 


1.47 






3 


3645 


0.52 


1.60 


1.71 


5.16 


2 


4400 


0.41 


1.03 


1.66 


4.33 


2 


4400 


0.27 


0.68 






6 


5500 


0.49 


1.03 


1.52 


3.19 


2 


7000 


1.14 


1.99 


1.86 


3.25 


2 


9000 






1.86 


2.66 


2 


12500 






1.77 


2.17 


1 


16200 






1.50 


1.74 


1 


22000 






0.65 


0.72 


1 


34000 






0.14 


0.14 


1 


erg cm" 


-^s-U- 


-1 









Menzies et al. 


|l982|) 3 


5 -ISchild et al. 


moi). 6 
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in which d is the distance to the object, V is the volume of the ionized zone 
and e\ (ergcm^ s~^ A~^) is the total (f-b + f-f) volume emission coefficient per 
electron and per ion in the wavelength scal e . We calculated them according to 
expressions given by iBrown and Mathewsl ()l97Cll ) adopting the Gaunt factors 
to be unity. In s pite of this approximation, our procedure re produces the 
tabulated values ( Brown and Mathewsl . 19701 : Gurzadvan . 1997 ) within ~5% 
at the blue side of the Balmer discontinuity and within 1% or less at other 
wavelengths for electron temperatures Te ~10 000 - 20 000 K. The difference 
is smaller for higher T^. 



3.2 Blackhody accretion disk and the boundary layer continuum 



Because of having defined the continuum only by fluxes at selected wavelengths 
(i.e. only the continuous radiation can be investigated), we will assume a 
steady-state, time-independent accretion disk around a central star of mass 
Mace and radius -Race- Further, the accr e tion d isk is optically thick and radiates 
locally like a blackbody (e.g. iTvlendal . Il98lh . The source of its energy is the 
gravitational potential energy of the accretor. The observed flux distribution 
of the disk. Fx, is then given by contributions of blackbody annuli integrated 
through the entire disk. The flux emitted by each annulus at the radial distance 
r is weighted by its area, 2nr dr. Thus the observed flux distribution of a disk 
at a distance d can be written as 



271 cosfil 



d^ 



Rd 

- J Bx{T,s{r))rdr, 



(2) 



where i is the inclination of the disk, its radius, and the radial temperature 
structure of the disk. 



Tes{r) = T, 



RsLCC 




Race 








- r ■ 




■ r ■ 



1/2-1 1/4 



(3) 



where 



T, = 4.10 X 10^i?9 ^/^(acc) Mi^^(acc) M^f (acc) K 



1/4, 



l>l/4^ 



(4) 



in which Rg a nd M^rr. are the radius and mass of the accretor and Mig is the 



accretion rate ( Warner . 1995[ l. The notation X„ means that the parameter X 
is expressed as fraction of 10" units (the units are cm, g,s). A useful character- 
istic of the disk is its maximum temperature, which occurs at r = (49/36)-Racc 
and has a value 0.488 T*. 
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Fig. 2. Parameter ^ as a function of the blackbody temperature. It represents a 
fraction of the hot star luminosity below 912 A (Sect. 4.2). 

The final stage of the mass transferred throughout the accretion disk is its 
deceleration and landing onto the accreting star. This happens in the so called 
boundary layer (BL), in which the inner disk material loses its kinetic energy 
to match the surface velocity of the star. If the accreting star rotates at a small 
fraction of its breakup velocity, then the BL emits a luminosity comparable to 
that of the disk. Contrary, when the star rota t es clo se to breakup, there should 
be little emission from the BL (e.g. Warner, 1995| ). The geometry of the BL 
can be approximated by a ring with thickness 2H surrounding the accreting 
star. Then the luminosity, Lbl, of its area, 27ri?acc2i?, for a non-rotating star 
is comparable to that of the accretion disk. So, 



BL 



1 GMaecM^^ 

2 Rficr 



(5) 



The effective ter nperatur e of th e BL can be derived with the help of Eq. 5, 
and according to IWarnei it can be approximated as 



Tbl ~ 2.9 x 10°M; 



1/3, 



acc) Rg '^^^(acc) Mig^l 



acc, 



K. 



(6) 



In the following section we use these components of radiation to fit the ob- 
served SED. 



4 Modeling the UV/optical continuum 



4-1 The uneclipsed nebula 



Figure 1 shows a faint component of the nebular radiation seen during eclipses. 
We compared a grid of models of the nebular continuum to these fluxes and 
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Fig. 3. Comparison of the observed ARPav continuum and calculated according 
to an ionization model: Hot star (T^ = 10^ K) ionizing the circumbinary medium 
giving rise to the nebular continuum (Te = 20 000 K). Solid thick line represents 
superposition of all components (G+H+NH-N2). These models are not able to give 
satisfactory agreement with observations (Sect. 4.2). 

selected that with Te = 40 000 ± 5 000 K, which corresponds to the best 
fit to observations (Fig. 1). Luminosity of this nebular emission is L-^ ~ 
500 ((i/4.9 kpc)^L0. The uncertainty in Te results from that in fluxes as given 
by the lUE observation s (about 10% , in average). An analysis of the x^{Te) 
function is described in Skopall ( 2001 ). 



4.2 SED of an ionization model 



During quiescent phases of symbiotic binaries the hot component (Th ~ 10^ K) 
ionizes a fraction of the neutral wind of the cool giant giving rise to the nebular 
continuum (Te ~ 1.5 — 2.0 x 10^ K). As a result the UV/optical continuum 
is given by superposition of the stellar and nebular components of radiation. 
Below we test this possibility for the case of ARPav. 

The maximum luminosity of a blackbody hot star of a radius i?h and temper- 
ature Th, which can be converted into the nebular radiation is 

912 

AtiRI J TiB^iT^) dX = CL^ erg s'S (7) 




where the factor (< 1) represents the fraction of the hot star luminosity, Lh, 
below 912 A. It is a strong function of the blackbody temperature between 
20 000 and 60 000 K (Fig. 2). Under assumption that the ionized medium is 
optically thick in the Lyman continuum, i.e. all the photons with A < 912 A 
are converted into the nebular emission, the maximum nebular luminosity can 
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be expressed as 

oo 

1"^:^= JexdxJn^n,dV = CL^ ergs'^ 



912 V 



Thus if the near-UV continuum is of the nebular nature, then for the ob- 
served flux, for example at A=3 600A, -F3^q^o ~ 5 x 10~^^ ergcm"^ s~^ 
(Table 1), the distance d = 4.9 kpc and the emission coefficient eseoo = 2.6 x 
10~^^ erg cm^ s~^ A^^ (Tc = 20 000 K), we get by using Eq. 1 the emission 
measure EM = Jyii^nedV = 5.5 x 10^°cm~^ and according to Eq. 8 

oo 

L^:^ = EM J exdX = l 230 Lq. (9) 

912 



The corresponding minimum of the hot star luminosity, which is capable of 
producing such nebular emission, is L™™ = L^^/( = 1380 for Th = 10^ K 
{( = 0.89). To compare the model with observations, first, we scaled the emis- 
sion coefficient e\ to fluxes between A3000-3600A, where the nebular com- 
ponent of radiation dominates the SED in the ionization model. The scaling 
factor /cn = F^^^/ex = EM/ind'^. Second, we scaled the hot star radiation 
with a factor > L™™/47r(i^ crTj^ to match observations in the far-UV. Figure 
3 shows our best solution, given by superposition of the nebular continuum 
with Te = 20 000 K and a hot stellar source (L^ ~ 2.5 x Lg^''^, Th = 10^ K), 
which roughly matches the observed UV continuum with a larger deviation at 
the far-UV. Cases with > 20 000 K signiflcantly exceed the far-UV fluxes, 
and those with Te <20 000 K cannot produce its rather flat proflle. All cases 
do not match the optical continuum. Finally, the hot star with a lower tem- 
perature transforms a smaller fraction of its stellar radiation into the nebular 
(e.g. Th = 510^K gives ( = 0.57), which then requires too high T™" other 
than to match the SED in the far-UV. 

This analysis shows that the observed SED of the AR Pav UV continuum does 
not correspond to that of the ionization model of symbiotic binaries during 
quiescent phases. On the oth er hand, the presence of highly ionized emission 
lines in the UV spectrum fiHutchings et al. . 1983[ l and the nebular emission 



during eclipses indicate a source of ionizing photons in the system. 



4-3 SEDs with accretion disk 



In accordance with previous suggestions, here we will flt the UV/optical/near- 
IR continuum of the hot object in ARPav by models with an accretion disk 
as introduced in Sect. 3.2. 
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4-3.1 SED with a single accretion disk 



Here we calculated a grid of models of fluxes given by Eq. 2 for fixed values 
of Mace and Mace, but different -Race, and selected that giving the best fit to 
the observed fluxes under conditions of the least square method. 



According to ISchild et al.l (|2001[ ) we adopted the mass of the donor (the giant 
star) Mg = 2M0 and that of t he accretor M^c = 0.75 Mtn. The rate of 



period change P = —3.5 x 10 ^ fjSkopal et al.L l2000a| ). the mass ratio q 



Mg/Macc = 2.7 and the orbital period of 605 days yield the accretion rate 
-^acc = 1.7 10~^ Mq ?/r~^ for the Roche lobe overflow if t he transferred mass 
and its angular momentum are conserved in the system ( Pringle and Wad3 . 



19851 see Sect. 5.1 for more details). Using these values of Mace and Mace to 



fit the component of radiation supposed to be due to the accretion disk (= 
SED - N - G; Fig. 4) by Eq. 2, we found the best solution for R^^ = 0.78 i?© 
(Fig. 4, Table 2). We integrated the disk contributions to its outer radius 
i?d = 10 Rq, which can be considered as the smallest dimension of the disk for 
these parameters. This is because the disk temperature between the surface 
of the accretor and r ~ 10 Rq is higher than 5 000 K, and thus radiative 
contributions of this part of the disk are important to explain the UV/optical 
SED (i.e. i?d cannot be smaller than 10 Rq), but the disk annuli with r > 10 Rq 
(i.e. T(r) < 5 000K) do not contribute significant emission at UV/optical 
wavelengths, and are negligible with respect to the radiation of the giant 
in the IR region (i.e. Rd can be larger than ~ 10 Rq). Superposition of the 
considered components of radiation here (G, N, AD) matches well the observed 
SED of ARPav (Fig. 4). This model expresses the round profile of the far-UV 
continuum, the fiat near- UV/optical continuum and a relatively significant 
(with respect to the ionization model) contribution to the near-IR continuum 
of the giant. Note that we observe about 0.2 mag deep eclipse even in the / 
band (|Menzies et al.l . ll982h . 



The accretion disk luminosity here. Lad ~ 1 700 ((i/4.9 kpc)^//©, can be bal- 
anced by the energy liberated in an accretion process for large accretion rates 
of Mace > 10~^ M0?/r~^ (Table 2). In addition, for the orbital inclination of 
~ 70° (Sect. 5.2), and assuming the disk plane to be coplanar with it, the 
total disk luminosity is by a factor of about 1.5 larger (L^d ~ 2 cos(z)Lad)- 
As discussed below (Sects. 5.1 and 5.2), observations do not contradict such 
a large Mace for a certain period of evolution of the giant star in ARPav. 



4-3.2 SED with the boundary layer and accretion disk 

To calculate the flux contributions from the BL we assumed that this region 
can be approximated by a blackbody at a constant temperature Tbl (Eq. 6). 
In our procedure we fitted superposition of fluxes from the accretion disk and 
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Fig. 4. Comparison of the ARPav continuum and the calculated flux distribution 
for an accretion disk (AD). Solid thick line represents the resulting continuum given 
by the sum of all components of radiation (G+N+AD). Dashed thick line as in Fig. 
1. This model is described in Sect. 4.3.1. 

the BL, 



F, = F,{AD)+7tB^{Tbl), 



(10) 



to the observed SED, where Fx{AD) is given by Eq. 2. First, we scaled con- 
tributions from the BL with a factor 

kBL^ jFx{AD)d\/aT^^ (11) 



to satisfy conditions of Eq. 5. Second, we calculated a grid of models for Mace 
= 0.75 and 1.0 Mq, = 30 i?©, M^c between 1.010"^ and 5.0 10"^ yr"^ 
and i?acc between 0.1 and 3.4 it!©- Finally, we selected models, for which the 
function 



reached a minimum. We used fluxes from 1270 to 7000 A (Table 1), but omitted 
those at 2200, 2600 and 2700 A, i.e. N = 11. The function x is plotted in Fig. 5. 
As can be seen from this figure there is no satisfactory solution. The best fits 
correspond to radii i?acc between 0.7 and 2.5 Rq and the BL temperature 
Tbl ~ 40000 - 60 000 K. The impossibility to obtain a better solution is 
reflected by the too low required temperature of the BL, of which the flux 
dominates the far-UV wavelengths. Figure 6 shows an example. The results of 
this section suggest that the BL emission is not important to fit the ARPav 
SED. Below, in Sect. 5.3, we discuss the BL emission in more detail. 
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Fig. 5. The quantity x &s given by Eq. 12 in modeling the observed SED by the BL 
and accretion disk contributions. Sohd hues represent models for Mace = 1 
and Mace = 0.1 — 5.1 x lO~^M0?/r^^. Compared are also two x functions for 
Maec = 0.75 M0 (dotted lines). The minima of x correspond to Tbl ~ 40 000 - 
60 000 K. Their large values reflect poor fits of these models to observations (Sect. 
4.3.2, Fig. 6). 
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Fig. 6. One example from the best solutions of fitting the SED by the accretion 
disk and its BL flux contributions. The fits are very poor, because of a significant 
contribution of the BL in the far-UV wavelengths. This model is described in Sect. 
4.3.2 



5 Discussion 



5.1 A mass transfer problem in ARPav 



The accretion rate of 1.7 10^^ yr^^ (Sect. 4.3.1) derived from the orbital 
period change corresponds to the conservative mass transfer from the giant 
star onto its companion. However, successful fits of the UV continuum require 
at least one order of magnitude higher rates to balance the observed luminosity. 
This discrepancy can be caus ed by mass loss from the system, which m akes the 
orbital period increase (cf. ISoberman. Phinnev and van den HeuveL . .1997.) . 
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Table 2 

Examples of models of the SED for selected values of Ma, 



-^acc 




Mace 


acc -^max / -^B L 


x/10-13 


[^0] 


[Re] 


[Mq 2/r-i] 


[^0] [K] 








Accretion 


disk (Sect. 4.3.1) 




3.7 


40 


1.710-3 


5410 16600/46100 


0.08 


1.7 


20 


1.710-^ 


1180 16000/50600 


0.09 


0.78 


10 


1.710-5 


257 16800/55600 


0.11 


0.01 


0.1 


3.310-11 


0.04 16500/88500 


0.07 


Accretion 


disk and the boundary layer (Sect. 


4.3.2) 


2.0 


30 


2.110"^ 


1240 15600/46700 


2.51 


2.5 


30 


3.610-^ 


1700 15100/44300 


2.53 



Notes: In aU cases Mace = 0.75 M©, Lace = l/2GMacc Macc/^acc- 

So, the mass transfer rate and the mass loss rate rival each other in the or- 
bital period change. Mass loss from ARPav is indic ated by r adial velocities 
of absorption components of hydrogen lines. Already Mavalll (1937) revealed 
a P-C ygni type of hydrogen line profiles on three plates from 1900, 1901 and 
1908. Thackerav and Hutching^ ( 1974 ) measured absorption components of 



the H7 and/or B.5 line profiles between 1953 and 1972. They determined an 



averag e position of these absorptions at -85.2 kms ^. Recentlv. iQuiroga et al 



f)2002l ) measured radial velocity of the central absorption component of Ha at 



96 kms 1 on spectra taken between 1988 and 2001. These valu es are shifted 
by -17 to -28 km s~i relativ e to the systemic velocity (-68kms~i. lSchild et al. 



200l[ iQuiroga et all , liool and thus refiect a mass loss from the system. Un- 



der such conditions the observed decrease in the orbital period requires a 
much higher mass transfer than that given by conservative case. Therefore, 
the mass transfer rate (i.e. the accretion rate) of 1.7 10~^ MQyr~^ can be 
considered only as a lower limit. This fact allow us to estimate the accretor 
radius -Race > 0.1 Rq for Mace = 0.75 — 1.0 Mq to balance the observed accre- 
tion luminosity irrespectively to the accretion disk model. Such the large i?aee 
precludes a white dwarf to be the accreting star in ARPav. 



5.2 Roche lobe overflow or a wind accretion? 



Observations and the model energy distributions discussed in this paper sug- 
gest that the sole energy source in ARPav is accretion. This view is also 
supported by the behaviour in its historical, 1889-2001, light curve, in which 
no large (5 - 7 mag) outbursts as in symbiotic novae have been observed (see 
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Skopal et all |2001a|). In addition, the light curve has characteristic features 
similar to those observed in CVs (|Mavall Il937l : ISkopal et all Hoooj, which 
suggests that the basic principles of mass transfer in ARPav are similar to 
those of CVs. These arguments suggest that the red giant in ARPav fills up 
(or is close to) its Roche lobe. Contrary, if the giant in ARPav underfills its 
Roche lobe, then the mass transfer is via accretion from a stellar wind. How- 
ever, this cannot be true. If this were the case, there would have been a huge 
(unreahstic) mass loss via the wind (at least of 10~^ Mgyr"^) to get the ac- 
cretion rate required to balance the observed luminosity of the hot object. In 
that case, the wind mass loss would make the orbital period increase, while 
the observations show that the period decreases. The latter can only be due 
to the Roche lobe overflow. 



Other observational results allow this possibility. For example, t he gia nt's 



radius inferred from eclipses, Rg/A = 0.30 ± 0.02 (|Skopal et al.L |2000a|), is 



equal to its volume-e quivalent Roche radius for the orbital inclination of 70° 
f Quiroga et al. . 2002[). The V-shape of t he eclipse profile supports this possi- 



bility (e.g. Fig. 8 of lSkopal et al.Ll2000a| ). Further, the 106-day periodic varia- 



tion with the amplitu de of a few x 0.01 mag observed in the visual light curve 
fjSkopal et al.L l2000a| ) could be understood as a response of the giant's radius 
to its mass loss. Generally, after losing an amount of mass, the giant will restore 
its hydrostatic equilibrium, which leads to variation in difference of the radius 
of the star and its Roche lobe. In the case of the presence of a deep convective 
envelope in the giant star, its radius rapidly increases and the star can run into 
a vio l ent mass loss instability (see Soberman. Phinnev and van den Heuvel . 



1997 : Shore. Livio and van den Heuvel Il994 for more detail). Such a tran- 
sient high mass transfer rate from the red giant to the hot component could 
be responsible for a highe r level of the AR Pav activity observed during 1900- 
01, 1935-36 and 1985-99 (|SkoDal et allEToOb.) . 



5.3 Boundary layer emission in ARPav 



The nebular component of radiation indicates the presence of a hot ionizing 
source in ARPav. In the models with accretion disk the ionizing source is 
represented by the BL. 

The nebular emission has a luminosity Ln > 500 {d/A.O kpc)^L0 (= uneclipsed 
nebula without emission lines), which places strong constraint on the BL lu- 
minosity and ionizing capacity. For parameters of the fit shown in Fig. 4 and 
described in Sect. 4.3.1, Lbl{= L^c) ~ 1200 and Tbl ~ 51 000 K (Table 
2). To balance the Ln luminosity, a fraction of the BL emission below 912 A 
has to be 

Lbl(A < 912 A) = C i^BL > 500 Lq. (13) 
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For Lbl = 1200 {d/A.9kpcfLQ, Eq. 13 yields the factor C > 0.42, which forces 
Tbl > 40 000 K (Fig. 2). Thus the BL characterized by these quantities can 
balance well the observed nebular luminosity, L^. 

However, our modeling the SED with the BL (Sect. 4.3.2, Table 2) does not 
provide satisfactory fit (Figs. 5, 6). This implies that we cannot see the hotter 
inner parts of the disk and that any radiation generated in the BL is absorbed 
and diffused in the disk in the direction to the observer, but it is free in 
directions to the poles to ionize there the circumstellar matter. This situation 
thus produces two components of radiation of a very different temperature 
regimes observed in the ultraviolet spectrum: A relatively cool continuum, 
which was investigated in this paper, and a hot ionizing source seen through 
the nebular emission. 



It is plausible that the lack of the BL emission in the continuum results only 
from a different real geometry of the accretion disk (probably more extended 
at its edge) than we adopted in our modeling (Eq. 2). So, if the disk is rather 
thick in its outer parts, then we can see essentially only its outer edge, which 
occults the BL at the disk's center for a high orbital inclination. Such a disk 
structure was suggested for the first time by IPlavec and Hubenvl (|l994^ for 
interacting binaries in which the accreting star in non-degenerate. 



5.4 Configuration of the hot object and the eclipse profile 



Our models of the observed SED (Sects. 4.1, 4.3.1) suggest that the hot object 
in ARPav consists of an accretion disk around a central star embedded in a 
hot corona. This basic configuration can be independently verified by the 
eclipse profile, which is a function of wavelength. In the following two sections 
we therefore explore its two main characteristics - the different eclipse depth 
and width in the UV/optical/near-IR region on the basis of our model with 
accretion disk. 



5.4-1 Depth of eclipses 



Already Hutchings et al. ( 19831 ) noted that there is no total eclipse below 
3600 A, and 70% of the hottest continuum (below 1400 A) is seen even at 
central eclipse. Recentlv ISchild et al. (|200l[) demonstrated this behaviour on 
the ultraviolet lUE light curves. 



Generally, the eclipse depth is given by the ratio of fluxes observed in and out 
of eclipse. In view of our models, it is determined by the ratio of the N+G to 
all considered components of radiation. In Fig. 7 we compared the modeled 
eclipse depth according to the solution shown in Fig. 4 and that observed in 



15 



o 



3.5 
3 

2.5 
2 

1.5 
1 

0.5 



-0.5 



AR Pav 


_ 








\ ■ 


U B V R I 



3.2 



3.4 



3.6 



3.8 



Log wavelength [A] 

Fig. 7. Wavelength-dependent depth of echpses measured in the dereddened 
UV/optical/near-IR hght curves (full boxes). The solid line represents calculated 
eclipse depth according to our model in Fig. 4. 
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Fig. 8. The ratio of the hot object fluxes, and Fqs, for the solution shown in 
Fig. 4. This characterizes the broadening of the eclipse profile (see text). 

the light curves (cf . IS child et al. . 2001 : Andrewsl . 1974 : Menzies et al. 



We can see that our model explains well the wavelength-dependent depth of 
eclipses. 



5.4-2 Width of eclipses 



Observations sho w a very broad eclip se profile in the far-UV {A^p > 0.2 — 0.3 of 
the orbital phase, Schild et al. . 200l[ ). while in the optical/near-IR the profile 
is characterized by a. narr ow core {A(p < 0.1) with shallow wings (l Andrewi . 



1974l:lMenzies et aLllTosj ). The hot object in ARPav consists of basically two 



components - the extended nebula and the central star with an accretion disk. 
Let /n and fcs be the fiuxes emitted by a surface element of the projected 
nebula and the central source, respectively, to the sky. The observed eclipse 
profile then can be characterized by the ratio /n//cs- The case of /n//cs ~ 1 
makes the eclipse profile very broad, because the extended nebula contributes 
significant emission with respect to that of the small central source. On the 
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other hand, /n//cs ^ 1 results in a narrow eclipse core determined by the 
central sourc e, accompanied with shallow wings as the nebula is faint (see 
also Fig. 4 of lSkopal et ID . bOOlbh . In our approximation and for a spherical 
nebula, it is simple to find that 



/n/ /cs oc Fn/ -Fcs, 



(14) 



where and Fcs are integrated fluxes of the nebula and the central source 
at a given wavelength. Figure 4 shows that the ratio F^/Fqs is a function 
of the wavelength, which thus represents the primary cause of the observed 
wavelength-dependent eclipse width. Therefore, we plotted this ratio in Fig. 8 
according to modeled components of radiation shown in Fig. 4. In the effect 
of relation (14), Fig. 8 shows that 

(/n//cs)uv > (/N//cs)opt, (15) 



which answers the question why the eclipses are broader in the ultraviolet than 
in the optical/near-IR region. A more rigorous approach, in which the param- 
eter /n//cs is treated as a function of the radial distance from the central 
object, would give a more accurate information about the nebula structure. 
However, this analysis is beyond the scope of this paper. 



6 Conclusions 



The results of this study may be summarized as follows: 

(i) Based on the ultraviolet lUE low-resolution spectra and the broad-band op- 
tical and infrared UBVRI and JHKL photometry, we re-constructed the SED 
in the ARPav continuum between 0.12 and 3.4 /im. By this way we revealed a 
very high luminosity of the hot object in the binary, Lh ~ 2 200 ((i/4.9 kpc)^ Lq. 

(ii) We failed to fit the UV/optical continuum of the hot object by a standard 
ionization model of symbiotic binaries. This suggests that the configuration of 
AR Pav differs significantly from the typical configuration of symbiotic binaries 
during their quiescent phases. 

(iii) The UV/optical/IR SED can be fitted by a simple blackbody accretion 
disk embedded in an extended hot corona (Te = 40 000 ±5 000 K). The stan- 
dard thin disk model requires a high accretion rate M^^c > 2 x 10~^ Mq yr~^ 
onto the central star with a radius -Race ^ 2Rq. Any accretion process at 
Mace > 1-7 X IQ-^M^yr-^ (Sect. 5.1) for Mace = 0.75 - 1.0 Mq limits the 
accretor radius to i?ace > 0.1 Rq to balance the observed accretion luminosity. 
This result suggests that the accreting star in ARPav is non-degenerate. 
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(iv) Te basic configuration of the hot object - the accretion disk embedded 
in an extended nebula - allow us to understand the observed wavelength- 
dependent depth and width of the eclipse profile (Figs. 7, 8): The eclipse 
width results from that the hot object consists of the two geometrically very 
different components of radiation (N + AD), the ratio of which is a function of 
the wavelength. The eclipse depth is given by a different ratio of the uneclipsed 
(N+G) to the total hght throughout the spectrum. 

(v) Models with the accretion disk and its BL (Lbl ~ -^ad) failed to fit the 
far-UV SED. However, the presence of a very hot region in the system is 
indirectly indicated through the strong nebular emission. This suggests that 
the hotter inner parts of the disk with the BL are permanently occulted by 
the disk material, when viewing the system (nearly) edge-on. 
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